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ABSTRACT 

We construct and analyze a u-band selected galaxy sample from the SDSS Southern 
Survey, which covers 275 deg 2 . The sample includes 43223 galaxies with spectroscopic 
redshifts in the range 0.005 < z < 0.3 and with 14.5 < u < 20.5. The S/N in 
the it-band Petrosian aperture is improved by coadding multiple epochs of imaging 
data and by including sky-subtraction corrections. Luminosity functions for the near- 
UV 01 u band (A rs 322 ± 26 nm) are determined in redshift slices of width 0.02, 
q • which show a highly significant evolution in M* of —0.8 ± 0.1 mag between z — 

and 0.3; with M* - 5logh 70 = -18.84 ± 0.05 (AB mag), log0* = -2.06 ± 0.03 
(h% Mpc" 3 ) and logp L = 19.11 ±0.02 {h m WHz" 1 Mpc" 3 ) at z = 0.1. The faint-end 
d ■ slope determined for z < 0.06 is given by a — —1.05 ± 0.08. This is in agreement with 

recent determinations from GALEX at shorter wavelengths. Comparing our z < 0.3 
luminosity density measurements with 0.2 < z < 1.2 from COMBO-17, we find that 
■ the 280-nm density evolves as pl oc (l + z) 13 with (3 — 2.1±0.2; and find no evidence for 

any change in slope over this redshift range. By comparing with other measurements 
of cosmic star formation history, we estimate that the effective dust attenuation at 
280 nm has increased by 0.8 ± 0.3 mag between z = and 1. 

Key words: surveys — galaxies: evolution — galaxies: fundamental parameters — 
galaxies: luminosity function, mass function — ultraviolet: galaxies. 



1 INTRODUCTION 

In the absence of dust, the rest-frame UV luminosity of a 
galaxy is nearly proportional to the total mass of short- 
lived OB stars and therefore to the star formation rate 
(SFR). This has been used to show that the volume- 
averaged SFR of the Universe h as been declining s i nce at 
least z ~ 1 JLilly et all Il996t : iMadau et alJ Il99rl Il99st 
ICowie et ailll999t ISteidel et al.lll999t IWilson et alJEoolh 
Until the observations of the Galaxy Evolution Explorer 
(GALEX) iMartin et al.ll2005ft . the accuracy in the mea- 
sured evolution rate had been low because of the lack of 
100-300 nm surveys at low redshift. The FOCA ballon-borne 



telescope iMilliard et al.ll 992f) had been the only instrument 
for measuring the galaxy luminosity function (LF) at these 
wave l engths; this survey co vered about 2 deg 2 llTrever et alJ 
1998; Sull ivan et alJ l2000f) . However, 'sun-tanning bands' 
u/U/U' also provide a window on the star-forming prop- 
erties of galaxies as they are significantly more sensitive 
to young stellar populations than the B/bj/g bands; over 
5000 deg 2 has been imaged in the it-band as part of the Sloan 
Digita l Sky Survey (SDSS). Using these data. lHopkins et alJ 
( 2003) have demonstrated that the u-band luminosity is a 
reasonable measure of the SFR by comparing with other 
SFR indicators (Ha, Oil and far-IR). 

Figure^shows the normalized responses of the GALEX, 
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Figure 1. Filter pr ofiles: GALE X iMartin et al.ll2003ft: FOCA 
<Milliard et alj|l992ft : SDSS uariz iStoughtorj l et alJl2002l) : APM 
fej; and an artificial 280-nm band JWolf et al.l2003T) . The u-band 
shifted to z = 0.1 is shown by the dash-and-dotted line ( 01 u 
band). Each curve is normalized so that the integral of the nor- 
malized transmission TdlnA is equal to unity. Thus the height 
of each curve is related to its resolving power (A/ A A). The thick 
gray line represents the difference in magnitudes between dust- 
free 5 My r and 1 0000 Myr si mple s tellar populations from the 
models of Bruzual & Chariot (2003) (the j/-axis for this line is 
shown on the right). 



FOCA and SDSS filters. The difference between the spec- 
tra of a young and an old stellar population is also plotted. 
This shows (i) the importance of the u-band in constrain- 
ing the spectral energy distributions of galaxies between the 
GALEX UV bands and the visible bands, and (ii) the sig- 
nificantly increased sensitivity to young stellar populations 
of the u-band in comparison with the g-band. 

While imaging data can be used to study correlations 
between properties of galaxies, in order to measure the 
space density of sources, it is preferable to select the tar- 
gets for spectroscopic followup using the most appropriate 
band. The largest spect roscopic surveys at z < 0.2 are the 
SDSS iYork et al.ll2000t). the 2 dF Galaxy Redshift Survey 
RS) (|C( 



. . , JColless et al.ll200lf) 
vey (6dFGS) Jjones et al.ll2004l) 



(2dFGRS) 



and the 6dF Galaxy Sur- 
which have principally se- 
lected galaxies using the r, bj and K bands, respectively 
(A e ft ~ 616, 456 and 2160 nm). While the bj selection is more 
biased toward star-forming galaxies than the other two, the 
band does not sample a majority of the light below the rest- 
frame 400- nm break at z < 0.15. In this paper, we describe 
the SDSS u-band Galaxy Survey (uGS), which provides an 
intermediate selection between the r/bj and the GALEX 
bands, and a local sample for comparison with higher red- 
shift surveys that are selected at simil ar rest-frame wave- 
lengths (e.g. DEEP2. iDavis et aljEoolh . The u-band inte- 
grates flux almost entirely from below the 400-nm break 
(A c ff ?s 355 nm). 

The plan of the paper is as follows: in §|5|we describe the 



The SDSS u-band has a small amount of contamination from 
around 710 nm. This corresponds to about a 0.02 mag effect fo r 
mid-K stars and galaxies of similar color (Abazaiian et al 12004) . 
We have not corrected for this because most of the galaxies in 
our sample are significantly bluer. 



basics of the SDSS (which can be skipped for those familiar 
with the main survey) ; in § [3] we introduce the Southern 
Survey; in §0]we present the results for the LFs (with more 
details in Appendix[XJ; in §§|K|&|S|we discuss and conclude; 
and in Appendix [B] we outline sky-subtraction corrections 
for SDSS u-band magnitudes. Note that magnitudes used 
in this paper are corrected for Milky- Way (MW) extinction 
unless otherwise noted. 



2 THE SLOAN DIGITAL SKY SURVEY 

The Sloan Digital Sky Survey is a project, with a dedicated 
2.5-m telescope, de signed to image 10 4 deg 2 and obtain spec- 
tra of 10 6 objects jYork et alj|2000l : IStoughton et alJl200l 
lAbazaiian et all2003U2004l) . The imaging covers five broad- 
bands, ugriz with effective wavelengths of 35 5, 467, 616, 
747 a nd 892 nm, using a mosaic CCD camera iGunn et alJ 
Il998j). Observatio ns with a 0.5-m photometric telescope 
jHogg et alll200ll) are used to calibrat e the 2.5-m telescope 
images using a st andard star system jFukugita et al.lll996l : 
ISmith et al]l2002l) . Spectra are obtained using a 640-fiber 
fed spectrograph with a waveleng th range of 380 to 92 nm 
and a resolution of A/AA ~ 1800 llUomoto et alJll999T> . 

The imagin g data are astrometrically calibrated 
iPier et alJ 120031) and the images are reduced using a 
pipeline photo that measures the observing conditions, and 
detects and measures objects. In particular, photo pro- 
duces various types of magnitude measurement including: 
(i) 'Petrosian magnitudes', the summed flux in an aper- 
ture that depends on the surface-brightness (SB) profile 
of the o bject, a modified version of the flux quantity de- 
fined by IPetrosianl lll976t) : (ii) 'm odel magnitudes ', a fit 
to the flux using the best fit of a de_Vaucouleursl <ll959l) 
and an exponential profile jFreeman ^7ol) : (hi) 'PSF mag- 
nitudes', a fit using the local point-spread function. The 
magnit udes are extinction-corrected using the MW dust 
maps of lSchlegel. Finkbeiner fc Davij (^)98j). D etails of the 
imaging pipelines are given by iLuoton^t alJ i200 ll) and 
IStoughton et alJ i2002l) . 

Once a sufficiently large area of sky has been imaged, 
the data are analyzed using 'targeting' software routines 
that determine the objects to be observed spectroscopically. 
The targets that are part of the 'main pro gram' include 
galaxies with r 



•p gr; 

pctro < 17.8 ('MAIN selection' Hstrauss et ail 
2002); quasars selected by vario us color criteria with i us { < 

19.1 or 20.2 ('QSO selection'; iRichards et"aD Eo02): and 
luminous-red galaxies selected by color cuts with r p ctr o < 

19.2 or 19.5 ('LRG selection': lEisenstein et~all 12001?) . The 
targets from all the samples are assi gned to plates, each wi th 
640 fibers, using a tiling algorithm jBlanton et aLll2003dl) . 



3 THE SOUTHERN SURVEY AND SAMPLE 
SELECTION 

The main program of the SDSS is concentrated in the North- 
ern Galactic Pole (NGP), with only three 'stripes' (2.5° 
wide) in the Southern Galactic Pole (SGP). 2 During the 

2 In the nomenclature of SDSS, a 'stripe' consists of a North- 
ern and a Southern 'strip' because there are gaps between the 
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times when it is not possible to observe the NGP, the 'South- 
ern Survey' has been in operation. This has involved re- 
peat imaging of the middle SGP stripe (on the celestial 
equator) and additional non-standard spectroscopic obser- 
vations. Here we define the Southern Survey region as RA 
from -50.8° to 58.6° (20.6 h to 3.9 h) and DEC from -1.26° 
to 1.26° (an area of 275.7 deg 2 ). Over this region, there are 6 
to 18 repeat images depending on the sky position with 253 
unique spectroscopic plates observed (including 57 for the 
main program); up to JD 24 53228. The SDSS uGS consists 
of high signal-to- noise ratio (S /N) magnitude measurements 
using a coadded catalog fS I3.1fl . sky-subtraction corrections 
to u-band Petrosian magnitudes (Appendix^J , and spectra 
selected in a variety of ways (S 13. 21 with completeness cor- 
rections for Up C t r o < 20.5 (§EHnj- These sections (§§02ES1 
Tabled Figs.|5J]SJ can be skipped or browsed for those not 
interested in the survey details at this stage. 



3.1 Creating a coadded imaging catalog 

The repeat imaging data can be coadded at the image level 
or at the catalog level. The former is necessary for increasing 
the depth of the imaging. Here we are mostly interested in 
improving the u-band Petrosian S/N and associated colors. 
Therefore, it is adequate to coadd the data at the catalog 
level because even the bluest galaxies in our sample are de- 
tected in the g and r bands with adequate S/N to define 
a Petrosian aperture. This has the advantage that the data 
passes through the standard imaging pipeline with no ad- 
justments. 

Our procedure for coadding the catalog was in two 
parts: (1) producing a master catalog of groups of matched 
objects; and (2) selecting appropriate samples for analysis. 

(1) First of all, each camera column and 'strip' is con- 
sidered separately (see footnote |5J • Catalog objects are se- 
lected from all the runs so that they are: (i) unique to a run 
('status_ok_run'); (ii) not near the edge of a frame (not 'ob- 
ject_edge') unless they are deblended from an edge object 
('object_deblended_at_edge'); (iii) not 'bright', which refers 
to a preliminary identification of bright objects in the cata- 
log; (iv) not 'blended' unless they are the final product of a 
deblending process ('object_nodeblend' or 'nchild'=0); and 
(v) detected in two or more of the the five bands, which 
is to remove cosmic rays and other artifacts only detected 
in one band. All the qualifying objects are then matched 
to objects in other runs within a radius of 1.5 arcsec. If 
there is more than one object from one run in a single 
matched group (which is rare but can occur because of de- 
blending), then only the object that has the smallest de- 
viation from the median r pc t ro value of the group is se- 
lected. Various mean and median values are calculated for 
the matched groups of objects, while the object flags are 
taken from the highest quality run. T he SDSS asinh magni- 
tudes llLupton. Gunn fc Szalavll999l) are converted to linear 



detectors, which are aligned in six camera columns. The two 
strips are interleaved to produce a contiguously imaged stripe. 
A ' run' is a continuous dr ift-scan observation of a single strip. 
See lStoughton et alJ i2002fl for details on nomenclature including 
the flags used to select the objects. 



units before averaging fEq. 11321 . This procedure produced a 
coadded catalog of 6 x 10 6 objects. 

(2) In addition to selecting on magnitudes, we apply the 
following criteria: (i) the object is detected in at least two 
bands out of ugr based on the best imaging run; (ii) the 
matched object is from a combination including over half 
of the maximum number of runs available at that sky po- 
sition; 3 (iii) the object is not saturated in the fiducial r- 
band; (iv) the S/N of the Petrosian flux is greater than three 
in each of the u, g and r bands (the uncertainty is deter- 
mined from the standard error over the coadded runs). The 
u-band Petrosian flux is adjusted for sky-subtraction errors 
(Appendix^ and if there are large differences between the 
median and mean ugr fluxes, the median is used to replace 
the mean value (and a modified standard error that excludes 
the minimum and maximum values is used). For the SDSS 
uGS, we selected resolved sources with ttpetro < 20.5. 

The data from the coadded imaging catalog that we use 
for the science results in this paper are: the mean Petrosian 
fluxes (median for ~ 3% of objects); the standard errors of 
the fluxes 4 (+2% error added in quadrature for fc-correction 
fitting); the median r ps f — r mo doi values (for star-galaxy sep- 
aration); the mean sky positions; and the mean Petrosian 
half-light radii. The median S /N of the coadded-catalog Pet- 
rosian u-band fluxes at u ~ 20.5 is 10, which is a factor of 
~ 3 improvement over a single epoch of imaging. 

3.2 Spectroscopic target selections 

In addition to repeat imaging along the Southern Survey 
equatorial stripe, there are also extra spectroscopic observa- 
tions. There are six general programs that contribute most 
of the redshifts to a u pct ro < 20.5 galaxy sample including 
specifically designed selection criteria for this survey. 

For the selection, star-galaxy separation and SB param- 
eters are defined as follows: 

Asg = ^psf fmodcl (1) 

A*r,50 = r po tro + 2.51og(27r7i 2 50 ) (2) 

where R r ,50 is the radius containing half the Petrosian flux. 
Thus u r ,5o is the mean SB within the Petrosian half-light 
radius. Note that most of the targets were selected using a 
single epoch of imaging before any coadded imaging catalog 
was produced. 

(1) MAIN selection. Many of the spectra were targeted 
as part of the main galaxy sample. This has the following 
basic criteria: 

r pctro < 17.8 (3) 
p r ,5o < 24.5 (4) 
A sg > 0.25. (5) 

3 Complex objects could have their centres shifted by more than 
the matching radius between runs and matching in more than 
half the runs ensures that only one copy of a complex object is 
included. 

4 The standard error is a/V~N; while the modified error is given 
by 1.3(T mo d/VN — 2 when using the median flux, where <r moc i 
is the standard deviation excluding the minimum and maximum 
values. The factor of 1.3 is determined from the median of 0"/o" mo d 
for the data. 
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Table 1. Number of objects and spectra for u-band selected samples 



sample selection 


no. of 


no. 


spectroscopically observed by SDSS (broken down by prof 


;ram) a 




2dF b 




objects 




T2)pSO 


pjlt-band 


|(4)|ow-z 


|(5)[ow-z 




other 


total 




u < 20.0 (incl. unresolved) c 


321768 


18829 


15061 


5606 


2888 


1027 


139 


8317 


51867 


79 


u < 20.5 & A sg > 0.05 d 


74901 


22141 


6825 


7400 


6687 


1846 


637 


303 


45839 


104 


u < 21.0 & A sg > 0.05 e 


146488 


23502 


10089 


7536 


9997 


2200 


2084 


584 


55992 


109 



"Number observed by various Southern Survey programs described in § 13. 21 and the total number for all the programs; up to JD 
24 53228. 

'Number of 2dFGRS redshifts not observed by SDSS. Only two 2dFGRS fields overlap with the Southern Survey. 
c Sample used to assess star-galaxy separation; see Fig. [3] for redshift versus A sg (Eg. HI. 

''Sample used to compute the galaxy LFs (§HJ; see Fig. |2]for the redshift histograms, Fig, l4lfor color-color plots, and Fig, IBIfor 
spectroscopic completeness as a function of color. 

e Sample used to assess completeness to fainter magnitudes; see Fig.|5]for completeness as a function of magnitude. 



See lStrauss et all l)2002|) for details. In addition to the main 
program, galaxy targets that were missed, because of fiber- 
placement restrictions or photometric errors, were included 
on additional plates. Thus the completeness of the MAIN 
selection is very high (« 99%) over much of the Southern 
Survey. 

(2) QSO selection. SDSS quasars were selected by looking 
for non-stellar colors using PSF magnitudes. For the low 
redshift candidates there was no requirement that the object 
be a point source. This enables selection of resolved galaxies 
where the central light may be dominated by a quasar but 
the total light may not be. The low redshift (i.e. extended 
sources) selection in the main program can be approximated 
by 



(4) Low-z galaxy selection. Targets were selected using 
photometric redshifts with primarily: 



2psf 

Upsf ,9psf 

S f — 0.5p ps f — 0.5r pst 



< 19.1 

< 0.9 

< 0.9 . 

See lRichards et alJ <2002l) for details (e.g. fig. 13 of 1 
per). In addition to the main program, several alternative 
target selections were made on additional plates. These in- 
cluded probing closer to the stellar locu s (and also closer to 
the g alaxy locus), in the stellar locus llVanden Berk et alJ 
2005) and selecting fainter QSO candidates. 

(3) u-band galaxy selection. Here the idea was to was to 
obtain a complete u-band magnitude-limited galaxy sample 
by filling in the redshifts missed by other selections. The 
criteria were mostly as follows: 

(9) 
(10) 

(11) 
(12) 
(13) 

where u sc icct = u mo dci-r modc i+r pct r , which can be regarded 
as a pseudo-Petrosian u-band magnitude. The reason for 
using this magnitude definition was to avoid u po t r o (single 
epoch), which has larger Poisson noise and systematic errors 
(Appendix [BJ , and has a greater susceptibility to imaging 
artifacts. The g and r-band requirements were also used to 
reduce artifacts. On some plates, the magnitude criteria were 
relaxed by 0.2 mag to u S eiect < 20.0. 



Uscicct 


< 


19.8 


^?pctro 


< 


20.5 


7"pctro 


< 


20.5 


Mr,50 


< 


24.5 


Asg 


> 


0.2 . 



Zphoto < 0.15 

rpctro < 19.0 

A SE > 0.15 . 



(14) 
(15) 
(16) 



The photometric technique was calibrated using spectro- 
scopically confirmed redshifts from the main program and 
the Southern Survey. In order to test selection effects, other 
photometric redshift ranges were also sparsely sampled as 
were galaxies with 19.0 < r pe tro < 19.5. 

(5) Low-z galaxy selection (precursor). Targets were se- 
lected using photometric redshifts with primarily: 



(6) 








(7) 


•2photo *C 


0.2 


(17) 


(8) 


Rpetro *C 


20.0 


(18) 


Asg > 


0.15 . 


(19) 


pa- 









•^photo 


> 


0.3 


^modcl 


< 


19.5 


Asg 


> 


0.15 



These plates were also designed to be complimentary to the 
higher photometric redshift selection described below. 

(6) High-z galaxy selection. Here the idea was to obtain 
spectroscopic redshifts for non-LRG galaxies above redshift 
0.3 in order to improve photometric redshifts for these types 
of galaxies. The targets were selected with: 



(20) 
(21) 
(22) 



In detail, the photometric redshift cut was converted to a 
series of color and magnitude cuts. 

Spectra were matched to photometric objects within 1.5 
arcsec or within R r .so (1.5-12 arcsec, which accounts for 44 
objects where an earlier version of photo may have targeted 
a deblended piece of the galaxy). The number of spectra 
contributing to various u-band selected samples are shown 
in Tabled Figure|5|shows the redshift histograms. Other se- 
lections include stellar programs for unresolved sources and 
about 100 redshifts were included from the 2dFGRS (which 
is similar to g < 19 galaxy selection). 
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Figure 2. Redshift histograms (ftpctro < 20.5 sample): number per 0.005 bin versus redshift. The top-left panel shows data from all the 
spectroscopic programs; with the solid line representing the observed da ta and the dotted line representing the completeness-corrected 
data 13.31 . The dashed line is a fit using Eq. 1 of Pcrcival ct al. (2001) with parameter values 0.158, 2.054 and 0.603. The other five 
panels show the data from the six programs described in § 13.21 Note the j/-axis scales vary; and the program names refer to the selection 
algorithms not the spectral classifications. 



3.3 Completeness corrections 

To determine spectroscopic completeness corrections, we do 
not attempt to back track and reproduce the above selection 
criteria. Many of the criteria are complicated and somewhat 
arbitrary for creating a u-band magnitude-limited sample. 
In addition, the photometric code and imaging run varies 
between selections and our aim is to select using Petrosian 
magnitudes from a coadded catalog. Instead we use an em- 
pirical approach for estimating the completeness factors (C), 
which we define as the fraction of photometric objects that 
have been observed spectroscopically. 

Before describing the completeness corrections, we dis- 
cuss star(-quasar)-galaxy separation. For each object we 
used the median A sg value (= r ps f — r mo( jei) from the coad- 
ded imaging for robustness. Figure [3] shows redshift versus 
A sg for a spectroscopic sample with M po tro < 20.0 (Tablc0. 
We define four regimes: sources with intra-galactic redshifts 
(z < 0.002); and unresolved (A sg < 0.05), weakly resolved 
(0.05 < A sg < 0.4) and strongly resolved (A sg > 0.4) 
sources with extra-galactic redshifts (z > 0.002). The main 
point to note is that there are no unresolved extra-galactic 
sources with redshifts less than about 0.2 in this sample. 
While there are certainly selection effects against targeting 
this type of object, we note that about 7000 of the point 
sources with redshifts were targeted at random (as part of 



testing the QSO selection: fVanden Berk et al.ll2005T) . There 
are some weakly resolved sources at low redshift. In this 
regime, some objects are spectrally classified as quasars 5 
and some as galaxies. 

To determine completeness corrections, we divide the 
sample into bins using four variables: A sg , u, u — g, and 
g — r. These are related to the primary selection variables de- 
scribed in S 13.21 First, galaxies are divided into two samples 
based on a cut in A sg at 0.3 (near the limit for MAIN selec- 
tion). Figure^Jshows u — g versus g — r distributions for the 
strongly resolved and weakly resolved sources (now divided 
at 0.3). The kinked dashed li nes show the limit fo r QSO 
selection of extended sources (Richards ct al 1 l2002l) while 
the straight dashed line shows a cut at U — r = 2.2, which 
approxi mately divides the b imodality in the galaxy distri- 
bution JStrateva et alJl200ll) . The point to note is that the 
completeness will vary with position in this diagram, be- 
cause of cuts in r, g or associated colors. In addition, photo- 
metric redshift selection will also depend strongly on these 
colors because the 400-nm break moves through the p-band 
over the redshift range 0.0 to 0.3. 



By spectral classification as a QSO, we mean a spectrum with 
broad emission lines, i.e., classified as a type 1 (unobscured) AGN. 
There is no consideration of luminosity or (narrow-)line ratios. 
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Figure 3. Redshift versus star-galaxy separation parameter for 51236 objects with u pe tro < 20.0. The dashed lines divide regions where 
the spectral classification is highly uniform (> 99.9% STAR; > 99.9% QSO; 99% GALAXY) except for the QSO+GALAXY region, 
which has a mixed classification (77% and 23%, respectively). The axes are linear in log(z + 0.004) and log(A sg + 0.1). The solid lines 
represent logarithmically-spaced density contours with 4 contours per factor of 10. Note that there are no point sources (A sg < 0.05) 
with redshifts between 0.002 and 0.195. There are two ID-spectral routines used for redshift and spectral classification within the SDSS 
collaboration (M. Subbarao et al.; D. Schlegel et al. in preparation). Only redshifts where the two routines produced similar results were 
included in this figure. This rejects only 1% of the measured redshifts. 



After dividing the sample using A sg , for u > 19.7, the 
strongly resolved objects are divided into 0.1-mag bins and 
the weakly resolved into 0.2-mag bins. At brighter magni- 
tudes, wider bins are used. These magnitude bins are further 
divided in u — g and g — r with from 2 x 2 to 16 x 8 color 
bins such that there are a minimum of fifty objects per final 
bin. Note that we have not included SB as a variable which 
would be necess ary for studying bivariate distributions in- 
volving SB (e.g. lCross et al-lboOltlBlanton et al.ll2003<J> . 

Figure shows the average completeness as function 
of magnitude for four groups divided by A sg and by color. 
The completeness is higher for the strongly resolved objects, 
with the redder objects having higher completeness because 
of MAIN selection. In the weakly resolved group, the bluer 
objects have higher completeness because of QSO selection. 
For the group dominated by late-type galaxies (solid line 
in Fig. EJ, the average completeness is 88% at u — 19.5 
dropping smoothly to 18% at u = 20.5. Above this limit, 
the completeness drops below 2% for some bins (A ag > 0.3). 
Figure HJ shows completeness as a function of color. 

For the SDSS uGS catalog, we use the following selec- 
tion criteria from the coadded imaging catalog (S I3.11 : 

Upetro < 20.5 (23) 



g pctI o < 21.0 (24) 
A sg > 0.05 (25) 

^extinction < 1.4. (26) 

The later cut on MW extinction only excludes objects over 
~ 0.7 deg 2 around RA of 58.4° . 6 We have not included an ex- 
plicit SB limit because the production of the coadded catalog 
(§ 13.11 is robust to low-SB artifacts. However, there is an im- 
plicit limit because of the pipeline reduction jBlanton et alJ 
120041) . This does not affect our results significantly, which 
are focused on the luminosity density (LD) measurements. 

The catalog includes 74901 objects, of which, 45839 
have been observed spectroscopically with SDSS (see TableQ 
for breakdown by program), and 104 with the 2dFGRS; the 
average completeness is 61%. The limit has been extended 
to 20.5 rather than stopping at the 19.8/20.0 limit of the 
M-band selection (item |(3)| of § 13.21 because pipeline flux 
changes and the other selections allow this. In other words, 



6 The area lost to bright stars is also small (< ldeg 2 in total), 
e.g., around the V = 3 star HR8414 and the globular cluster 
M2. We assume the area of the survey is 275 deg 2 for calculating 
galaxy number and luminosity densities. 
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Figure 4. Observed color-color plots for (a) strongly resolved sources and (b) weakly resolved sources. The kinked dashed lines represent 
Eqs. ITIandlSl Objects to the left of these could be selected by standard QSO selection (10% and 76% of those observed are spectrally 
classified as QSO, for each panel, respectively). The straight dashed line [in Panel (a)] represents an approximate division of the bimodality 
in the galaxy distribution (u — r = 2.2). Galaxies to the right are generally early types and those to the left are generally late types. 
The solid lines represent logarithmically-spaced density contours with 4 contours per factor of 10. Note that, for this plot, objects were 
restricted to those with S/N > 5 in u, g and r, and with ^extinction < 1- 




o.o L 

15 16 17 18 19 20 21 

Upetro 



Figure 5. Average completeness versus u-band magnitude for different types of objects. The solid and dashed lines represent strongly 
resolved objects that are bluer and redder than u — r = 2.2, respectively (Fig.|4ji); the dash-and-dotted and dotted lines represent weakly 
resolved objects that are bluer and redder than the approximate QSO cut, respectively (Eqs, ITIandlSl Fig. HbV 




completeness percentage 




Figure 6. Completeness as a function of g — r versus u — g for (a-c) strongly resolved sources with various magnitude ranges and (d) 
weakly resolved sources. Note there are no areas with < 1% completeness so white areas represent < 5 objects per 0.05 X 0.05 color bin 
[< 2 objects in Panel (d)] . The dashed lines represent the same color cuts as in Fig. |I] which shows the number densities. 
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all galaxy spectral types are sufficiently sampled spectro- 
scopically as faint as u ~ 20.5. Note that even at this limit, 
the redshift reliability is very high (~ 99%) because galax- 
ies are either bright enough in the visible (r pct ro < 19) or 
they very likely have strong emission lines. We assume that 
all the measured redshifts (D. Schlegel et al. pipeline) are 
correct for calculating the LFs. 



4 GALAXY LUMINOSITY FUNCTIONS 

From the wGS catalog fEqs. galaxies are selected 

with spectroscopic redshifts in the range 0.005 < z < 0.3 
and with magnitudes in the range 14.5 < u < 20.5. In ad- 
dition, we remove weakly resolved sources (A ag < 0.3) that 
are spectrally classified as a QSO (see footnote |SJ. In other 
words, only galaxies where the integrated visible flux is dom- 
inated by stellar light are included. These selections produce 
a sample of 43223 galaxies. Note that compact non-QSO 
galaxies (A 3g < 0.3) contribu t e < 1% to the LD. 

Following iBlanton et all (l2003rl . we fc-correct to the 
rest-frame band equivalent to the observed M-band at 
z = 0.1, which is called the 0A u band (A e ff ~ 322 nm; 
FW HM » 53 nm) . The absolute magnitude on the AB sys- 
tem llOke fc Gunnlll983ft is given by 

M 322 = itpotro - fco.it.,,. - 5 log(Di,/10 pc) - 0.04 (27) 

where: Dl is the luminosity distance for a cosmology with 
(fi m ,n A )o = (0.3,0.7) andtfo = (h 70 ) 70 km s ~ 1 Mpc~ 1 , and 
kg.iu.u is the k-cor rection using the m ethod of Blanton et alj 
l|2003al) (see, e.g.. lHogg et al.ll2002l for a general definition 
of the fc-correction) . The —0.04 term is the estimated cor- 
rection from the SPS S u-band system to an AB system 
jAbazaiian et al.ll2004li . 7 

To calculate the galaxy LFs, we divide the sample into 
0.02-redshift slices and 0.2- magnitude bins. For each bin, the 
LF is then given by 

<^dM = V — j (28) 

i 

where: C is the spectroscopic completeness (S I3.3L and V max 
is the c omoving volum e over which the galaxy could be ob- 
served <ISchmiddll96ct) within the redshift slice and within 
the magnitude limits of 14.5 < u < 20.5. Each redshift slice 
is cut in absolute magnitude so that the sample is nearly 
volume limited (0.9 < Vmax/V B u cc ^ 1 for ~ 95% of the 
galaxies) . The volumes of each slice range from 2 x 10 4 Mpc 3 
(z = 0.005-0.02) to 8 x 10 6 Mpc 3 (2 = 0.28-0.30). 

Figure [3 a) shows the binned galaxy LFs for the 15 
redshift slices out to z = 0.3. The dominant effect is that the 
function shifts toward higher luminosities with increasing 
redshift. 

The general form of the LFs can be parameterized using 
the lSchechterl £l976) function, which is given by 

(l)LdL=r ^y e -^^ (29) 

7 The SDSS photometry was originally designed to be calibrated 
to an AB scale but because of filter variations between the natural 
system of the telescop e and the standard star system it has been 
modified lAbazai ian et ai1l2003t) . 



where: 4>l&L is the comoving number density of galaxies 
with luminosity between L and L+dL; L* is the 'characteris- 
tic luminosity' (for the exponential cutoff); (f>* is the 'charac- 
teristic number density'; and a is the 'faint-end slope'. This 
equation is converted to magnitude form, with M* as the 
'characteristic magnitude', and the slope is then — 0.4(a+ 1) 
in a log number density versus magnitude plot. We use this 
parametric form to quantify the evolution in terms of M* 
and <f>* . 

Before quantifying the evolution, we determine the 
faint-end slope. The density of galaxies fainter than —17 
is only determined for the three lowest redshift slices. On 
the assumption that the faint-end slope does not vary, we 
fit the best-fit slope over these lowest slices (2 < 0.06). The 
Schechter fits are shown in Fig. Efb). The best-fit faint-end 
slope is given by 

a = -1.05 ± 0.03 [±0.07] (30) 

where: the first uncertainty is the standard error from the 
Poisson noise; and the second uncertainty is an estimate of 
the systematic error, obtained by comparing with the best 
fits using combinations of three out of the four lowest red- 
shift slices (Fig. IAH . These variations may reflect the fact 
that the Schechter function is not a perfect match to the 
LFs (and therefore depends on the magnitude range of the 
fitting), there may be systematic uncertainties because of 
photometric deblending or astrophysical changes because of 
large-scale structure. 

Assuming that the faint-end slope does not vary signifi- 
cantly with redshift, we can then characterize with more pre- 
cision the evolution of the exponential cutoff shown clearly 
in Fig. |7| To do this, we restrict the fitting to a between 
— 1.15 and —0.95. For each redshift slice, M* and (j>* are 
fitted marginalizing over the allowed range of a (Table lATt. 
For the four highest redshift slices (2 > 0.22), \og(j>* is con- 
strained to be greater than —2.2; and for the lowest redshift 
(2 < 0.02), M* is constrained to be brighter than —18.4. 
This is because these parameters are not realistically fit in 
these regimes. 

Figure |H| shows the evolution in the Schechter parame- 
ters. There is a highly significant detection of evolution in 
M * with a slope of —3.1 ± 0.2 per redshift. This means that 
the exponential cutoff becomes more luminous by about 0.9 
mag between redshifts 0.0 and 0.3. There is a marginal de- 
tection of evolution in (j>* , which is more strongly affected 
by cosmic variance. This is only a 1.5-cr detection using the 
redshift slices up to 0.22. The fits to the evolution are given 
in Table |21 Note that the results (for M* in particular) de- 
pend significantly on assumptions about a, for example, for 
a = —1.05 strictly fixed we obtain an M* evolution slope 
of -2.5 ±0.1. Thus, M* gets more luminous by 0.8 ±0.1 
mag (2 = 0.0 to 0.3) depending on the details of the LF 
evolution. 

The lowest redshift slice was not included in the M* 
evolution fit. This was to avoid problems with measur- 
ing bright nearby galaxies using the SDSS pipelines. Large 
galaxies may be too strongly deblended (SB fluctuations are 
more significant at low redshift), or the Petrosian aperture 
may not be large enough. In these cases, the flux of the 
galaxies will be underestimated. 

The Schechter function allows for an estimate of the 
total comoving LD assuming that the function remains valid 
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absolute magnitude: M322 - 5 log h?o absolute magnitude: M322 — 5 log b.70 

Figure 7. Near-UV °- 1 u LFs for galaxies with 0.005 < z < 0.3. The lines represent different redshift slices. The binning is 0.02 in redshift 
and 0.2 in magnitude. There is a clear increase in the comoving density of luminous galaxies with increasing redshift. Panel (a) shows 
the binned functions, while Panel (b) shows the binned-function error bars and Schechter fits for the three lowest redshift slices and two 
others (fixed a). The full sample over the redshift range includes 43223 galaxies, which is reduced to 39494 by making the slices nearly 
volume limited. 




Figure 8. Evolution in the u LFs characterized by the Schechter parameters M* and 0*. The faint-end slope a is assumed to be in 
the range —1.15 to —0.95. The squares with solid error bars represent standard errors while the dotted error bars [in Panel (b)] represent 
standard errors plus an error added in quadrature to account for cosmic variance (Appendix^}. The triangles with dashed lines represent 
the 1-cr ranges when the fitting if restricted to M* < —18.4 or log0* > —2.2. The thick solid lines represent fits to the evolution over 
the range shown (Table l2l . The dotted line [in Panel (a)] represents a fit with fixed a = —1.05. 
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Table 2. Straight line fits to evolution in the Schechter parameters and the luminosity density 



parameter (322-nm band) redshift line intercepts line slope /3 value" 

fitting range (z = 0.0) (z = 0.1) (2 = 0.2) (per unit z) 

M* — 5 log ^70 (mag) b 0.02-0.30 -18.53±0.02 -18.84 ±0.01 -19.15±0.02 -3.06 ±0.18 3.3±0.2 

log(</>* /Mpc" 3 ) - 31ogh 70 0.00-0.22 - 2.02 ± 0.04 -2.06±0.02 -2.11 ± 0.03 -0.46 ±0.31 -1.2 ±0.8 

j + 2.51og/t 70 (mag Mpc -3 ) ... 0.00-0.30 -13.47±0.07 -13.68±0.03 -13.89±0.06 -2.12 ±0.55 2.2 ± 0.6 

M* — 5 log ^70 (mag) c 0.04-0.30 -18.64±0.02 -18.89 ±0.01 -19.14 ±0.01 -2.51 ±0.12 2.7 ±0.1 

\og(<t>* /Mpc" 3 ) - 31ogh 7 o 0.00-0.22 -2.05 ± 0.04 -2.08±0.02 -2.11 ±0.02 -0.27±0.29 -0.7 ±0.8 

j + 2.51og?t70 (mag Mpc -3 ) ... 0.00-0.30 -13.53±0.06 -13.71 ±0.03 -13.89± 0.03 -1.80 ±0.37 1.9 ± 0.4 



a Fit to linear measurements (L* , <f>* and pi,), with a function oc (1 ± z)@ . 

6 For the first set of fits, the results were marginalized over the faint-end slope: a in the range —1.15 to —0.95. By comparison with fixed 
a results, estimates of the systematic uncertainties are 0.05 in M* , 0.02 in log tj>* and 0.04 in j for the intercepts. 

c For the second set of fits, the faint-end slope was fixed: a = —1.05. Note here the 0.02—0.04 slice is an outlier in M* (Table |A*T1 . which 
may reflect photometric errors or large-scale structure variations. 
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Figure 9. Evolution in the ' u comoving LD. See Fig. GO for 

symbol and line-style details. 



outside the magnitude range of the fitting. This LD is given 
in magnitudes per cubic megaparsec by 

j = M* - 2.5 log [OVMpc -3 ) r f (a + 2)] (31) 

where Tf is the gamma function; and in linear units by 



ph = 10 O4.1-i)/2.5 WHz -l Mpc 



(32) 



from j in AB mag Mpc -3 . If we assume that 4>* an d a 
are constant then the evolution in M* also represents the 
evolution in the comoving LD. However, we cannot rule out 
contributions from variations in <j>* (or a) to the evolution. 
Figure [5] shows the evolution in the LD that includes the 
variation in <j>* . Parameterizing the evolution as 

p L <x(l + zf (33) 

as per for example iLillv et al] lll996j) . then we obtain f) = 
2.2±0.6 for the 322-nm LD. In § 15.21 we combine these mea- 
surem ents with COMBO-17 at 0.2 < z < 1.2 iWolf et al] 
2003). 



5 DISCUSSION 

5.1 Comparison with z ~ 0.1 luminosity functions 

The firs t test we mak e is to compare with the SDSS re- 
sults of iBlanton et al] (l2003bl) that used the same band. 
For the 01 u LF, their sample included 22020 galaxies with 
u < 18.4 in the redshift range 0.02 < z < 0.14. The u 
limit was chosen in order to avoid significant bias from the 
r < 17.8 MAIN selection, i.e., because 99% of galaxies have 
u ~ r ^ 0.6 [Fig. Ufa)]. The fitting used a maximum likeli- 
hood method, with a general shape for the LF, that incorpo- 
rated luminosity and number evolution. The results obtained 
(by fitting a Schechter function) were M* — —18.70 ± 0.03, 
a = -0.92 ± 0.07 and j = -13.71 ± 0.14 at z = 0.1 (after 
converting to Ho = 70). The faint-end slope and LD are 
in 2-<t statistical agreement with our results while M* is 
not (0.14 mag difference). However, if we set a to be —0.92 
in our analysis, the discrepancy is reduced. In addition, we 
used fixed 0.2-magnitude bins whereas they used multiple- 
Gaussian fitting. Note also that there are less than 4000 
galaxies in common between the two data sets (< 20% of 
theirs and < 10% of ours). 

Parameters were used to include evoluti on in the lu- 
minosity and number density ]Lin et aljll999l). In terms of 
the evolution fitting given in Tablel2l lBlanton et alf s results 
can be considered as —4.2 ± 0.9 for M* and +1.3 ± 1.3 for 
\og(j>* (both per unit z). These evolution parameters are in 
agreement with our results within 2-a. 

Using the uGS, we have improved on the accuracy of the 
01 u LF and evolution compared to using the main galaxy 
sample selected in the r band (which was of course a strong 
motivation for this survey). For example, the 01 u LD is now 
known t o comparable accurac y as the 01 griz bands mea- 
sured bv lBlanton et al] ]2003bl) . Figure tTol shows the z = 0.1 
luminosity densities from the SDSS. The uncertainties are of 
order 5-10% because of absolute calibration, SB selection ef- 
fects, conversions to total galaxy magnitudes and estimates 
of the effective survey areas. 

Figure 1101 also shows the l uminosity densities from 
the ballon-borne FOCA telescope JSullivan et al.ll200d) and 
from the GALEX space telescope jBudavari et al]l2005l) . 
The former was derived from 2.2 deg 2 over 0.0 < z < 0.4 
and the latter from 44 deg 2 over 0.07 < z < 0.13. The ear- 
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Figure 10. Luminosity densities at 2 = 0.1 from SDSS, FOCA 
and GALEX. The square represents the 01 u result of this pa - 
per; the triangles, the 01 griz r esults o f [Blanton^^all l2003hTl : 
the cross, the FOCA result of ISullivan et al .1 <2000l): and the 
diamonds, the GALEX results of iBudavari et al . |2005|. The 
ho rizontal bars represent the FWHM of each band. The best- 
fit iFioc fc Rocca-Volmerangd il99l Il999l) pop ulat ion-synthesis 
model from the fitting of TBaldrvfcGiazebrookl J2003fl is shown 
by the dotted line; and a new best-fit model that uses the GALEX 
results is shown by the dashed line. 



lier result is significantly higher in LD. Also. ISullivan et alJ 
obtained a faint-end slope of —1.5 ± .1, which is not in 
agreement with our measurements, while lBudavari et alJ ob- 
tained —1.1 ± 0.1, which is in good agreement. In Fig. 1101 
we also show two spectral models f or the UV to near-IR 
' cosm ic spectrum' from the fitting of Baldry & Glazcbrook 
(2003) (using PEGASE models: IFioc fc Rocca-Volmerangg 
ll997l.ll999T) . We use the new fit to the GALEX plus SDSS 
measurements to obtain a correction from the 0A u band to 
a 280 ± 20 nm band (Fig.0, which we use to compare with 
higher redshift measurements of rest-frame UV densities. 
This is given by 



J280 « J322 + 0.45 



(34) 



A similar correc tion is determined if w e use the fc-correction 
templates from iBlanton et al] (l2003aTl . i.e., fc-correcting to 
~ °' 3 u band, but this is less reliable because the z ~ 0.1 LF 
then requires significant wavelength extrapolation from the 
observed bands. 

Our low-redshift °u results (Table IAH are also in 
good agreement with that determ i ned in nearby clusters by 
IChristlein. Mcintosh fc ZabludofiU2004l) . They found My = 
-18.9 ± 0.3 (converting to AB mag) and au = -IT ± 0.2. 



5.2 



Evolution in comoving luminosity densities 

z ~ 1, we take 



near-UV LD measurements from the literature 


Lillv et al. 


1996 
2003 


; JConnollv^^^J 
; Budavari ct al. 


1997; Wilson et al] 2002; 


Wolf et al. 


2005s: IWvder et all 


2005) and convert 



them to the (Cl m , Q.a)o = (0.3,0.7) cosmology. This was done 
by making a comoving volume correction over the redshift 



range of each measurement and a luminosity correction for 
the midpoint redshift. Figure FiTI shows these measurements 
and our results converted to 280 nm, which show an increase 
in the LD with redshift. Our results and COMBO-17's re- 
sults (both shown using solid-line error bars) cover a sig- 
nificantly larger volume than the other results (excepting 
the GALEX results). We fit to these SDSS and COMBO-17 
results, and obtain 



/3 2 so = 2.07 ±0.14 [±0.10] 



(35) 



where the first uncertainty is the standard error and the sec- 
ond is per 0.05 mag uncertainty in the Eq. 1341 conversion. 
This represents the most accurate measurement of near-UV 



found bv iLillv et al. 


( 1996J) and is closer to t 


rates determined by 


Cowie. Soneaila & Barecr 



and 

IWilson et all <2002l) . The Eq.[|3fit is shown by a solid line 
in Fig. Illf b). while a fit using pl oc exp(t/r), where t is 
the look-back time for our assumed cosmology, is shown 
by a dashed line. The exponential timescale is given by 
T280 = 5.5 ± 0.4 [±0.3] Gyr. 

In the absence of dust and chemical evolution, the 280- 
nm density evolution corresponds closely to the SFR density. 
From a solar-metallicity population synthesis model (with 
a f3 ~ 2 star- formation history): /3sfr ~ /?280 ±0.1 and 
tsfr ~ i"280 — 0.2 Gyr. In other words, the evolution in the 
SFR is slightly steeper than that of the near-UV LD because 
of relatively small contributions from evolved stellar popu- 
lations. O ur results are consistent wi th that found using the 
far-UV bv lSchiminovich et al] ]2005l) : /3 150 = 2.5 ± 0.7. 

Many other SFR indicators have been used to trace the 
cosmic star formation history, some of which are less sensi- 
tive to dust or can be corrected for dust. From a compila- 
tion of UV, [Oi l], Ha, H/3, mid -IR, sub-mm, radio and X-ray 
measurements, Hooking (2004) found that /3sfr = 3.3 ± 0.3 
for z < 1 data when including SFR-dependent dust attenu- 
ation corrections where necessary. This is inconsistent with 
our result /3sfr = 2.2 ± 0.2 that assumes no dust evolution. 
The results can be reconciled if the effective average atten- 
uation at 280 nm increases by 0.8 ± 0.3 mag between z — 
and 1. This is a plausible increase because the characteristic 
luminosities o f rapidly star-form ing galaxies are increasing 
with redshift dCowie et aljll996l) ; and more l uminous galax- 
ies h ave higher levels of dust attenuation dHopkins et al] 
l200ll) . However, we caution that our measured dust in- 
crease is not independent of t hat idea because of the SFR- 
dependent corrections used bv lHopkinsl Bolometric OB-star 
luminosity densities derived from far/near- UV and mid/far- 
IR wavelengths (with corrections for AGN and evolved stel- 
lar populations) could be used for an unambiguous measure 
of Psfr- 

5.3 Evolution in luminosity functions 

While it is robust to compare our LD results with those of 
COMBO-17 using a model cosmic-spectrum correction, it is 
not so straightforward for assessing LF evolution. Any color- 
magnitude relations and dispersion will affect the number 
density and shape of the LFs. Instead, the COMBO-17 data 
could be analyzed to measure the LFs in the 0,1 tt or the more 
standard °u band (Appendix [XJ . Other large surveys for 
which our results could be compared are the VIMOS-VLT 
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Figure 11. Evolution in the near-UV comoving luminosity densities. Various data measurements are shown by the symbols. Panel (a) 
shows the data at wavelengths 230, 250, 280 and 322 nm; while Panel (b) shows the data ad justed to 280 nm where necessary. The 
corrections in logp^ were assumed to be —0.18, +0.18 and +0.14 for the 01 u, GALEX NUV and lWilson et all results, respectively (from 
the dashed-line fit in Fig. 1101 . Thus, the agreement between the results from this paper and from GALEX in Panel (b) is by construction. 
The best fit to the data from this paper and from COMBO-17 (fig. 19 of Wolf et all) is shown by a solid line using the parameterization 
of Eq . 1331 and by a dashed line using an exponentially increasing LD with look-back time. The x-axis is linear in log(l + z). 



Deep S urvey (VVDS: iLe Fevre et al]|20ol) or the DEEP2 
survey iDavis et al.ll20ojh ' 

Compa r ing ou r °u results with the VVDS results of 
lllbert et alJ i2004f) . we find M* brightens by 1.5 to 2 mag 
between z — and 1. Figure WI\ shows the evolution in M* 
versus redshift. To obtain upper and lower limits, we fit to 
the SDSS data (±2<r) and the VVDS data (±lcr) allowing 
for the errors to all be in the same direction. There is a sug- 
gestion that M* brightens more rapidly below z = 0.5 than 
above but the data are also consistent with a constant slope. 



6 CONCLUSIONS 

We have analyzed a it-band selected galaxy survey, us- 
ing spectroscopic completeness corrections as a function of 
color and magnitude to account for inhomogeneous selection 
(Figs. |5| |K|and |SJ. The main results are: 

(1) Testing star-galaxy separation (Fig. [3J, we find that 
compact galaxies that could be missed by MAIN selection 
contribute insignificantly to the stellar UV LD (< 1%). 

(2) The faint-end slope of the low-redshift 01 it LF is near 
flat: a = —1.05 ± 0.08. This was obtained from the best fit 
over the three lowest redshift slices (z < 0.06) with magni- 
tudes from -21.4 to -14.4 [Fig. Eft))]. 

(3) The evolution in the LFs is dominated by a luminosity 
shift (Figs. El and |HJ, which can be characterized by a shift 
in M* of -0.8 ±0.1 mag between z = and 0.3. 

(4) In order to compare the UV LD evoluti on with 
COMBO-17 at higher redshifts dWolf et alJ 12003ft . we fit 




redshift 

Figure 12. Evolution in M* for the °u LFs. The squares and 
diamonds represent the results from this paper, marginalized over 
a = —1.15 ± 0.10 and fixed a = —1.15, respectively (with errors 
approximat ely the size of sym bols). The triangles represent the 
results from lllbert et af] J2004I) with a allowed to vary except for 
the highest redshift slice. The sloped-line region shows the range 
between the lower and upper straight line fits. 



to LD measurements at z = 0.1 (Fi g. 1 1 lit including re- 
cent GALEX results of lBudavari et ail (120051) . and obtain a 
correction to 280 nm of +0.45 mag. Our LD measurements 
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versus redshift then lineup remarkably well with those of 
C0MB0-17 (Fig. [TTJ and we find that the evolution can 
be parameterized by Ea. I33l with /feso = 2.1 ± 0.2. This is a 
significantly shallower evolution than that found using other 
SFR indicators and is consistent with an increase in average 
dust attenuation of 0.8 ± 0.3 mag between z — and 1. 

Future work with the uGS could include: assessing the UV 
LF evolution for different classes of objects (e.g., AGN, 
late/early- type galaxies), measuring the Ha LF (modulo 
aperture and dust corrections) and studying multi-variate 
distributions in a star-forming sample. These will enable a 
quantification of contributions to the near-UV decline at 
z < 0.2 from different galaxy populations and from dust 
and star formation. 
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APPENDIX A: MORE DETAILS ON THE 
LUMINOSITY FUNCTIONS 

We use a straightforward U max in slices approach for de- 
termining the LFs (§ EJ. Schechter functions are fitted to 
the LFs using standard least-squares routines. The errors 
on each 0.2-mag bin are taken to be a modified Poisson er- 
ror where the variance without weighting would be N + 2. 
This is an appropriate variance for low number counts (when 
the expected value is not known). Figure lA"Tl shows the best- 
fit faint-end slopes versus redshift. As there are significant 
degeneracies between M* and a, we considered a limited 
range in a, —1.15 to —0.95, in order to assess the evolution 
m M* and </>*. Table lA"Tl gives the Schechter parameters for 
the redshift slices along with some details of the LF fitting. 




0.00 0.05 0.10 0.15 0.20 

redshift 



Figure Al. Best-fit faint-end slopes for the 01 u LFs. The hori- 
zontal bars represent the redshift ranges while the vertical error 
bars represent the 1-a errors. The sloped-line region shows the al- 
lowed a range when fitting the other Schechter parameters. Note 
that we do not interpret the change in a from z = 0.0 to 0.1 as 
being caused by galaxy evolution. The change could be caused by 
the varying magnitude limits in the fitting (Table IA1I . For the 
lowest redshift slices there is a higher weight from luminosities 
fainter than ~ M* + 1; whereas at z ~ 0.05—0.1 the fitting is 
giving higher weight to the 'knee' of the LF. In other words, the 
Schechter function is not a perfect match to the LF as there ap- 
pears to be a change in slope around M322 ~ —16.5 (Fig.[7J. For 
the purposes of this paper, assessing the evolution in M* and j, 
we use the three-parameter Schechter function and compromise 
on a. Note also there could be systematic errors at low redshift 
because of deblending issues that increase the number of faint 
galaxies at the expense of bright galaxies. 



The top half of the table gives the absolute magni- 
tude ranges used in the fitting, the number of galaxies, the 
mean of the correction factors and the Schechter parame- 
ters for a marginalized over the range —1.15 to —0.95. The 
column l//j, bs gives the 1-a range (not including cosmic 
variance) of the inverse of the fraction of LD that is within 
the magnitude limits, i.e., it represents the extrapolation 
factor from the observed to total LD. The <j)* an d j errors 
include an estimate of cosmic variance, that is proportional 
to V -0 ' 3 , added in quadrature to the Poisson errors. This 
approximate power and the normali zation were estimated 
from fig. 3 of ISomerville et al.l (l2004h . The cosmic variance 
on j is assumed to be less than that on 0* because the spe- 
cific SFR per galaxy increases in low-density environments, 
which partly offsets changes in number density. 

The lower half of Table lATl gives the Schechter param- 
eters for fixed a = —1.05 both for the 01 u and the °u band. 
From the M* and j differences between the different bands, 
we obtain 

M355 ~ M322 - 0.35 (Al) 

for the average galaxy, which is in good agreement with that 
determined from the cosmic spectrum fit (Fig. 1101 . 

We provide these Schechter parameters to make com- 
parisons with other surveys but they are strictly only valid 
from about M* — 2.5 (or the bright limit) to the faint limit 
of the fitting. At magnitudes brighter than about M* — 2.5, 
there is a significant excess of galaxies above the Schechter 
fit (from a composite LF with each redshift slice shifted by 
M*). 
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Table Al. Luminosity function fitting and Schechter parameters a 



range of 


M322 limits 


no. of 


mean of 


mean of 




M322 paras. 


marginalized 


over a = - 


-1.05 ±0.10 


z 


bright 


faint 


galaxies 


I/O 


^slice/^max 




—M* 


- log 0* 


-3 




1/ fj, obs 


0.005-0.02 


-18 


.6 


-14.4 


601 


1.07 


1.01 


< 18.58 6 


2.05 ±0.15 


13.44 ± 0. 


.28 


1.32-1.55 


0.02-0.04 


—20 


.0 


-15.4 


2421 


1.17 


1.01 


18 


.61 ± 0.09 


2.10 ±0.09 


13.46 ± 0. 


.15 


1.06-1.14 


0.04-0.06 


— 21 


. 1 


-16.4 


4738 


1.15 


1.01 


18 


.64 ± 0.04 


1.96 ±0.06 


13.73 ±0. 


.11 


1.09-1.13 


0.06-0.08 


-21 


.4 


-17.0 


5334 


1.23 


1.01 


18 


.75 ± 0.02 


2.03 ±0.05 


13.64 ±0. 


.09 


1.17-1.21 


0.08-0.10 


-21 


.6 


-17.6 


5058 


1.25 


1.01 


18 


.84 ± 0.06 


2.12 ±0.05 


13.55 ±0. 


.09 


1.33-1.42 


0.10-0.12 


-21 


.6 


-18.0 


5143 


1.34 


1.03 


18 


.88 ± 0.02 


2.10 ±0.04 


13.60 ±0. 


.07 


1.46-1.51 


0.12-0.14 


-21 


.8 


-18.4 


5549 


1.47 


1.04 


18 


.90 ± 0.06 


1 99 ± 04 


13 95 + 


.09 


1.81-2.05 


0.14-0.16 


-22 


.0 


-18.8 


3549 


1.51 


1.04 


19 


.03 ± 0.06 


2.15 ±0.04 


13.74 ±0. 


.11 


2.25-2.64 


0.16-0.18 


-22 


.2 


-19.2 


2719 


1.49 


1.06 


19 


.06 ± 0.07 


2.11 ± 0.04 


13.82 ±0. 


,14 


2.95-3.74 


0.18-0.20 


-22 


.2 


-19.6 


1987 


1.48 


1.03 


19 


.13 ±0.07 


2.07 ± 0.05 


13.99 ±0. 


.17 


4.35-5.85 


0.20-0.22 


-22 


.2 


-19.8 


1259 


1.74 


1.06 


19 


.17 ±0.06 


2.12 ±0.05 


13.92 ±0. 


.19 


5.46-7.62 


0.22-0.24 


-22 


.4 


-20.2 


467 


1.71 


1.04 


19 


.16 ± 0.08 


> 2.02 c 


13.90 ±0. 


.29 


11.7-19.7 


0.24-0.26 


-22 


.4 


-20.4 


313 


1.92 


1.07 


19 


.27 ± 0.06 


> 2.09 c 


13.92 ±0. 


.22 


16.0-24.0 


0.26-0.28 


-22 


.4 


-20.6 


220 


2.17 


1.04 


19 


.36 ± 0.06 


> 2.10 c 


13.99 ±0. 


.21 


21.4-31.4 


0.28-0.30 


-22 


.4 


-20.8 


136 


2.17 


1.18 


19 


.44 ± 0.07 


> 2.05 c 


14.12 ±0. 


.27 


29.6-47.9 



range of M322 paras, with fixed a = —1.05 M355 paras, with fixed a = —1.05 d 







— A/* 




- log 4>* 




-3 




1/ fj, ohs 


—M* 




- log <t>* 




-3 








,obs 


0.005 


-0.02 


< 18.45 


b 


1.96 ±0.15 


13 


.53 ± 0. 


27 


1.50- 


-1.64 


< 18.86 


6 


2.00 ±0.15 


13. 


85 ± 


0. 


27 


1.44- 


-1.58 


0.02- 


-0.04 


18.47 ± 0. 


.05 


2.02 ±0.08 


13 


.44 ± 0. 


15 


1.04- 


-1.11 


18.88 ±0. 


.05 


2.06 ± 0.08 


13. 


76 ± 


0. 


15 


1.04- 


-1.11 


0.04- 


-0.06 


18.76 ±0. 


.03 


2.02 ±0.06 


13 


.75 ±0. 


11 


1.11- 


-1.15 


19.16 ±0. 


.03 


2.04 ±0.06 


11. 


08 ± 


0. 


11 


1.08- 


-1.11 


0.06- 


-0.08 


18.85 ±0. 


.02 


2.08 ±0.05 


13 


.67 ±0. 


09 


1.20- 


-1.24 


19.25 ±0. 


.02 


2.12 ±0.05 


13. 


99 ± 


0. 


.09 


1.20- 


1.24 


0.08- 


-0.10 


18.87 ± 0. 


.02 


2.14 ± 0.05 


13 


.57 ±0. 


08 


1.38- 


-1.42 


19.26 ±0. 


.02 


2.16 ±0.05 


13. 


88 ± 


0. 


08 


1.39- 


-1.43 


0.10- 


-0.12 


18.95 ±0. 


.02 


2.13 ±0.04 


13 


.65 ± 0. 


07 


1.54- 


-1.58 


19.32 ±0. 


.02 


2.15 ±0.04 


13. 


98 ± 


0. 


07 


1.56- 


1.60 


0.12- 


-0.14 


18.91 ± 0. 


.02 


2.00 ±0.04 


13 


.96 ± 0. 


07 


1.91- 


-1.98 


19.28 ±0. 


.02 


2.01 ± 0.04 


14. 


28 ± 


0. 


.07 


1.96- 


-2.02 


0.14-0.16 


19.01 ± 0. 


.02 


2.13 ±0.04 


13 


.71 ±0. 


06 


2.32- 


-2.43 


19.40 ± 0. 


.02 


2.17 ± 0.04 


14. 


01 ± 


0. 


06 


2.36- 


-2.48 


0.16- 


-0.18 


19.06 ±0. 


.03 


2.11 ± 0.04 


13 


.81 ± 0. 


07 


3.19- 


-3.41 


19.40 ± 0. 


.03 


2.10 ±0.04 


14. 


18 ± 


0. 


07 


4.30- 


4.70 


0.18- 


-0.20 


19.13 ±0. 


.03 


2.07 ± 0.04 


13 


.98 ±0. 


08 


4.76- 


-5.27 


19.43 ± 0. 


.04 


2.02 ±0.05 


14. 


.42 ± 


0. 


09 


7.57- 


8.73 



"See text of Appendix lAl for details. The units of M* , 4>* and j are AB mag, Mpc — 3 and AB mag Mpc - 3 , respectively. 

6 The M* limits shown for the lowest redshift slice are 1-cr limits when M* is constrained to be brighter than —18.4 and —18.8 for 

322 nm and 355 nm, respectively. 

c The rp* limits for the redshift slices above 0.22 are 1-cr limits when log</>* is constrained to be greater than —2.2. 
d The faint M355 limits for the fitting were (—0.2, —0.4, —0.6), with respect to the M322 limits, for the redshift ranges (0—0.06, 
0.06-0.16, 0.16—0.20), respectively. The bright M355 limits were —0.2 or —0.4 with respect to the M322 limits. Beyond z = 0.2, the °u 
LFs are significantly less reliable because of V m ax and fc-corrections, and a g-based sample would be more appropriate. 



APPENDIX B: SKY-SUBTRACTION 
CORRECTIONS FOR SDSS J7-BAND 
PETROSIAN FLUXES 

The equivalent depths across the SDSS magnitudes for 
galaxy target selection correspond to (u, g,r,i, z) ~ (19.8, 
18.6, 17.8, 17.4, 17.1). In other words, there are an equiv- 
alent number of galaxies per square degree (~ 100 deg -2 ) 
for galaxies brighter than these limits in each of the 
bands. These also approximately correspond to the aver- 
age spectral energy distribution of galaxies near those mag- 
nitudes. Comparing these limits (plus average MW extinc- 
tion) to the point-sour ce 95% completeness limits given in 
IStouehton et all i2002T) . we obtain differences of (1.9, 3.4, 
4.2, 3.8, 3.3). Thus, the u-band galaxy measurements are, on 
average, the closest to the survey imaging detection limit. 
In addition, systematic errors can dominate because scat- 
tered light can be significant in comparison with the sky 
flux (which is less of a problem in the z-band, for exam- 
ple). This can affect flat fielding and the zero point of the 



sky-subtracted frame and therefore is most significant for 
galaxies with low surface brightnesses. The distribution of 
u-band surface brightnesses in the Petrosian aperture (2x 
the Petrosian radius in the r-band) is shown in Fig. IBll for 
a galaxy sample. 

The variation of the average u—g galaxy color, as a func- 
tion of camera column, is shown in Fig. IB2l for different SB 
bins. There is a significant non-astrophysical variation of the 
average color and the amplitude of the variation increases 
with decreasing SB. The variation is reduced by using model 
colors because the effective aperture is smaller. The peak- 
to-peak systematic variation of u — g color is around 0.3-0.4 
for Petrosian magnitudes (lower two panels of Fig. IB2t and 
around 0.1-0.2 for model magnitudes, for galaxies with u- 
band surface brightnesses between 27 and 28 mag arcsec -2 . 

On the assumption that the dominate systematic error 
is due to sky-subtraction errors, we determined linear offsets 
in SB as a function of pixel position that minimized the 
color variation. This was done by using the three lowest SB 
bins and fitting a cubic polynomial to the implied offsets 
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Figure Bl. Distribution of u-band surface brightnesses for a 
galaxy sample with r pc t ro < 19.4. The vertical dotted lines show 
the positions of the cuts used to divide the sample for testing. 
Note that the surface brightnesses are for the SDSS Petrosian 
aperture, which is two times the Petrosian radius in the r-band, 
and the values have not been corrected for MW extinction (unlike 
all other flux measurements in this paper). 



Table Bl. Polynomial coefficients for the u-band Petrosian flux 
corrections as a function of objc_colc 



camcol 


Po/10- 12 


Pi/10" 15 


P2 / lO" 18 


Ps/10" 21 


1 


0.187 


-2.976 


3.767 


-1.101 


2 


0.902 


5.103 


-6.910 


2.063 


3 


-0.616 


-1.674 


2.391 


-0.749 


1 


0.485 


-1.027 


1.460 


-0.456 


5 


0.481 


-6.971 


7.650 


-2.181 


6 


0.242 


-1.933 


2.514 


-0.907 



See Eos. IBlHB2l for how to apply the correction, and see Fig. IB3 
for an illustration of the offsets for each of the camera columns. 
Note that these corrections were determined for photo v. 5.4 
and a Southern Survey coadded catalog. Similar offsets are 
observed using earlier versions of the pipeline and for different 
regions of the sky. 



for each camera column separately. The M pc tro fluxes were 
redetermined and the process was iterated until the results 
converged. The mean flux offset was normalized to zero. The 
polynomial coefficients are given in Table IHIl and should be 
applied so that 



Wflux = ""flux + P(objc_Colc) TV (2R r 



(Bl) 



where: ttfl ux is the pipeline Petrosian flux in units of 'mag- 
gies' [— 2.5 log(maggies) = mag]; P is the polynomial func- 
tion of the pixel position, which is different for each 'camcol'; 
and R r is the Petrosian radius in the r -band. Note that th e 
SDSS databases use asinh magnitudes llLupton et al.l [l999). 
which can be converted to flux by 



Xfl u 



sinh 



ln(10) 
2.5 



X n 



In(b) 



21, 



where the b values are (1.4, 0.9, 1.2 1.8, 7.4) 
ugriz, respectively dStouehton et al1l2002l) . 



xl(T 



(B2) 



for 



A plot of the offsets is shown in Fig. IB3l and the varia- 
tion in the u — g color after correction is shown in Fig. IB4I 
After correction, there is no significant variation with cam- 
era column for each SB bin (reduced % 2 ~ !)■ This implies 
that our assumption that sky subtraction was the dominate 
systematic error varying with camera column was valid. For 
galaxies with low surface brightnesses (28 mag arcsec -2 over 
the aperture), the correction can be up to ±0.3 mag depend- 
ing on the pixel position. 
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Figure B2. Variation in the average observed u — g galaxy colors as a function of the pixel position ('objc_colc') and the camera column 
('camcol'). The panels represent different SB bins. The reduced \ 2 values on the assumption of no variation with camera position are 
shown in the panels, with the dotted lines showing the mean color. Note that the SDSS cameras operate in drift-scan mode so that any 
systematic error only depends on the detector column (and not on the detector row). 
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Figure B3. Estimated flux offsets required to correct sky-subtraction errors in the u-band. The lines represent the polynomial fits to the 
data given in Table||n] The units are linear in flux where 10 12 maggies corresponds to a magnitude of 30 [i.e. — 2.5 log(maggies) = mag] . 
The average of the flux offsets was normalized to zero. Note that these offsets were determined using a Southern Survey coadded catalog 
where there was some smoothing over objc.colc (FWHM~50-200) because some drift-scan runs were offset slightly in DEC. 
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Figure B4. Same as Fig. IB2| cxccpt the u-band fluxes have been corrected using the functions given in Table |B*T1 
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